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Abstract. It is widely accepted that individual Galactic globular clusters harbor two 
coeval generations of stars, the first one bom with the 'standard' a-enhanced metal 
mixture observed in field Halo objects, the second one characterized by an anticor- 



O , related CN-ONa abundance pattern overimposed on the first generation, a-enhanced 

^^ ' metal mixture. We have investigated with appropriate stellar population synthesis mod- 



els how this second generation of stars affects the integrated spectrum of a typical 
metal rich Galactic globular cluster, like 47 Tuc. Our main conclusions are: 1) the age- 
sensitive Balmer line, Fe line and the [MgFe] indices widely used to determine age, 
Fe and total metallicity of extragalactic systems are largely insensitive to the second 
'.^ , generation population; 2) enhanced He in second generation stars affects the Balmer 

f^ ' line indices of the integrated spectra, through the change of the turn off temperature 

CO , and the horizontal branch morphology of the underlying isochrones, which translate 

into a bias towards slightly younger ages. 
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1. Introduction 

X 

H 

C^ ■ A large body of spectroscopic data published in the last 10 years has conclusively established 

the existence of primordial surface chemical abundance variations of C, N, O, Na in individual 

Galactic globular clusters (GCs) of all metallicities, from very metal rich objects like NGC6441, 

to the most metal poor ones like M15 (see, e.g.. Carretta et al. 2010') . 

The accepted working scenario to explain these CNONa anticorrelations prescribes that the 
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stars currently evolving in a GC were born with the observed CNONa patterns. Intermediate-mass 
asymptotic giant branch (AGB) stars in the range ~ 3 - SM© (and/or the slightly more massive 
super- AGB stars) and massi ve rotating stars are consider ed viable sources of the necessary heavy- 



element pollution (see, e.g. JVentura & D'Antonal2005h . Slow winds from the envelopes of AGB 



(or super-AGB stars) or from the equatorial disk formed around fast rotating massive stars inject 
matter into the intra-GC medium after a time of order 10^ - 10^ yr, depending on the polluter 
Provided that a significant fraction of the material is not lost from the cluster, new stars (second 
generation, but essentially coeval with the polluters' progenitors, given their short evolutionary 
timescales) may be able to form directly out of pristine gas polluted to varying degrees by these 
ejecta, that will show the observed anticorrelation patterns. Another by product of this pollution 
may possibly be an e nhanced initial He-abundance for these second generation stars. Recent 
numerical models by ID'Ercole et al.l (120101) have started to explore in quantitative details this 



broad picture. 

2. Models 



We have started by modelling a reference set of a-enhanced SSPs - metal mixture with [a[Ps\~ 
0.4 - with ages t=12 and 14 Gyr, which we identify hereafter as standard models. We have con- 
sidered in this investigation a single [Fe/H]=-0.7 (Z=0.008), typical of the metal rich subpopu- 
lations of Milky Way - and of 47Tuc, whose integrated spectrum is a benchmark for population 
synthesis models - and an initial He mass fraction 7=0.256. For both ages we have then con- 
sidered a second generation population {modified CNONa models) whose metal composition has 
C decreased by 0.30 dex, N increased by 1.20 dex, O decreased by 0.45 dex and Na increased by 
0.60 dex with respect to the first generation a-enhanced mixture, all other metal abundances being 
unchanged. This pattern is typi cal of valu es close to the upper end of the observed anticorrelation 
patterns in Galactic GCs (Carr etta et al.l l2009). The metal distribution of this second generation 
coeval population has the same C-(-N-i-0 sum and the same Fe abundance (as a consequence also 
the total metallicity Z will be practically the same) as the first generation composition, in agree- 
ment with spectroscopic measurements on second generation stars within individual Galactic GC. 
For both the 12 and 14 Gyr second generation populations we have accounted for two alternat- 
ive values of Y, e.g. 7=0.256 - as in the first stellar generation - and 7=0.300, to include a 
possible enhancement of He in second generation stars. This latter, hereafter identified as He en- 
hanced models, i s consistent with const raints on the typical enhancement of He in Galactic GCs, 
as determined bv lBragaglia et al.l (l2010l) . 



We stress that these choices for the chemical composition of second generation populations 
are very general, and should give us a realistic estimate of their impact on GC analyses us- 
ing SSPs. Overall, our selection of representative extreme values of the CNONa variations will 
provide a first important indication of the maximum effect of these abundance anomalies on GC 
integrated spectra, and will serve as a guideline to interpret the abundance pattern derived from 
fitting SSPs to the observed spectra of Galactic and extragalactic GCs. 



The SSP models were computed adopting isochrones from BaSTI jPietrinferni et al .1120061) 
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Figure 1. The spectral stellar population models computed for this investigation are shown in blue, red 
and brown lines (standard, modified CNONa and He-enhanced modes, respectively - see text for details). 
Models shown are those for 12Gyr. In gray line in shown the observed integrated spectrum of 47Tuc 
from ISchiavon et al.1 (l2005h . The model spectra were convolved to match the spectral resolution of the 
observation, FWHM ~ 3.1 A. 



and a stella r synthetic grid computed with ATLAS 12 (iKurucz & Avretllll98ll:ICastellill2005h and 



SYNTHE ( Kurucz & Avrettll98illSbordone et al.l2004 codes. One synthetic grid was computed 
for each set of models - reference a-enhanced, second gereration with Y=0.256 and second 
generation with Y=0.300 - properly taking the respective chemical mixture into account. From 
the appropriate isochrone and grid of synthetic stellar spectra we have cal culated the inte grated 
spectrum for these first and second generation populations, employing the iKroupal (1200 ih initial 
mass function. The integrated spectra modelled cover the wav elength region 3500 to 6000 A at a 
spectral resolution R - A/SA = 10000. We refer the reader to lCoelho. Percival & Salarisl(l201 lb 
for more details on our models. 



We show in Fig. 1 the spectral models for 12 Gyr populations. 



3. Results 



In Figs. 2 to 4 we zoom some regions of the spectra - known to be sensitive to C, N, Na and 
atomic O - to highlight the effect (or absence of an effect) of the abundance change on the stellar 
population spectra. In the case of C and N (Fig. 2) and Na (Fig. 4), the effect is clearly seen in 
the spectra. In the case of atomic oxygen, the effect is so small (lower than 0.5% in flux) that the 
two models are indistinguishable. 
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Figure 2. Standard and modified CNONa models are shown in blue and red respectively, for a spectral 
region sensitive to molecular carbon and nitrogen. Observed 47Tuc spectrum is shown in dark gray. 
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Figure 3. Same as in Fig. 2, for a spectral region sensitive to atomic oxygen. The ratio between the standard 
and modified CNONa model is shown in the bottom panel. 
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Figure 4. Same as in Fig. 2, for a spectral region sensitive to sodium. 



The effect of the CNONa modifi ed abundances on Lick/IDS spectral indices has been presen- 
ted in detail in ICoelho et al.l (1201 ih . To summarize, the presence of a second generation of GC 
stars with unchanged He-content affects appreciably only the Ca4227, G4300, CNi, CN2 and 
NaD metal indices. Very importantly, the variation of Ca4227 goes in the direction to mimic a 
lower Ca abundance, if this index is used as a measure of the Ca content. Our results support 



Lee & Worthevi (12005 ) suggestion that the effect of a second generation population with CNONa 
anticorrelations may explain the discrepancies they find when comparing Ca4227, CNi, CN2 and 
NaD index strength (on the Lick IDS system) from their a-enhanced SSP models, with data for 
Galactic and M31 GCs. On the other hand, age, Fe-abundance and Z inferred from Hy6-Fe5406 
and Hy6-[MgFe] diagrams (or the H(5f and, to a slightly smaller extent, the Hyp counterpart) are 
confirmed once again to be robust and insensitive to the chemical abundance pattern of second 
generation stars. 



When we consider the case of a second generation population with enhanced He-abundance 
(y=0.300), the values of all metal indices in this second generation population are essentially 
identical to the 7=0.256 case, but the Balmer line indices are changed by the increase of He, 
mimicking an age younger by 2 Gyr We detected that the changes in H/?, H(Jf and Hyp compared 
to the 7=0.256 case are due to differences in the underlying isochrone representative of the He- 
enhanced second generation population. The TO and main sequence of the underlying isochrone 
are hotter by ~100 K compared to the 7=0.256 isochrone. Given that the TO mass at fixed age 
is smaller (by about 0.07 M© in our case) in the He-enhanced isochrone, and we assume that the 
total mass lost along the RGB is the same (0. 1 1 M©), the mass evolving along the HB is smaller, 
hence the typical Tefj of HB stars is higher, again by ~ 100 K. These higher Tej values for TO and 
HB increase the value of the Balmer line indices. 



It remains yet to be studied what is the effect of modified CNONa (a nd He) abundances on 
the m etalhcites and ages derived via spectral fitting codes (e.g .Koleva. Prugniel. Bouchard & Wu 
20091) . On the other hand, one might test the possibility to measure abundances - beyond iron 
and a-elements - with spectral fitting (e.g. by using these models to differentially correct stellar 
populations, see lWalcher. Coelho. Gallazzi & Charlot.2009l) . 
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Figure 5. Spectral ratio between the modified CNONa and standard models. The strong features in the blue 
are related to CN bands. 
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Figure 6. Spectral ratio between the He enhanced and modified CNONa models. The change of the con- 
tinuum slope and the features related to H lines are a result of the hotter TO temperatures of the He enhanced 
model. 



It is likely possible to detect the existence of those abundance changes in the residuals of a 
spectral fitting analysis. At a first approximation, we could assume the residuals due a significant 
change on the CNONa abundance will show a pattern similar to the differential spectrum between 
our standard and modified CNONa models, shown in Fig. 5. The signal of C and N are clearly 
seen, as the modified CNONa abundances enhances the CN features by up to 30% in flux. In 
fact, primordial abundance variations such as the ones explored in this modelling might be an 
explanation for the long term problem of the particularly st rong CN bands in globular clusters 
in M31, when com pared to the ones of the Milky Way (e.g. IBurstein. Faber. Gaskell & Krumm 
1984tlTragerll2004 '). In Fig. 6 is shown the ratio between the He enhanced and modified CNONa 
models, to illustrate the effect of enhancing He alone. We note that features related to the hydro- 
gen lines and the small change of the continuum slope, characteristic of an integrated "hotter" 
model, are in agreement with the explanation above on the efifect of the He enhancement on the 
TO and HB temperatures. 
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4. Conclusions 

We explored the effects that primordial abundance changes of C, N, O and Na have on integrated 
spectra of a typical metal-rich galactic globular cluster These variations affect appreciably only 
the Ca4227, G4300, CNi, CN2 and NaD metal indices. Ages, Fe-abundances and Z inferred 
from Hy6-Fe5406 and Hy6-[MgFe] diagrams (or the H^p and, to a slightly smaller extent, the Hyp 
counterpart) are confirmed once again to be robust. If the CNONa changes are accompanied by 
He enhancement, there might be a bia s towards younger age determinations of up to 2 Gyr We 



refer the reader to lCoelho et al.l(l201 11) for the complete results 



The effect of these abundance changes on ages and metallicities derived via spectral fitting 
techniques remains to be quantified. It might be particularly interesting to explore if these prim- 
ordial abundance variations can explain the strong CN features observed in globular clusters in 
M3 1 (the author thanks Scott Trager for pointing this out). 
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